First-Principles Determination of the Structure of Magnesium Borohydride 
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The energy landscape of Mg(BH4)2 under pressure is explored by ah initio evolutionary calculations. 
Two new tetragonal structures, with space groups PA and IAi/acd are predicted to be lower in 
enthalpy, by 15.4 kJ/mol and 21.2 kJ/mol, respectively, than the earlier proposed PA^nm phase. 
We have simulated X-ray diffraction (XRD) spectra, lattice dynamics, and equations of state (EOS) 
of these phases. The density, volume contraction, bulk modulus, and the simulated XRD patterns 
of I4i/acd and P4 structures are in excellent agreement with the experimental results. 



PACS numbers: 61.50.Ks, 61.05.cp, 63.20.D- 

Light weight metal borohydrides have recently received 
much attention owing to their high gravimetric and 
volumetric hydrogen densities compared to other com- 
plex hydrides [Iffi]. Of these, magnesium borohydride, 
Mg(BH4)2, is a prominent lightweight solid-state hydro- 
gen storage material with a theoretical hydrogen capacity 
of 14.8 wt %. Hence it has been extensively studied as 
a template for developing novel hydrogen-storage solu- 
tions. Based on the experimental data, the ground state 
a and f3 phases, have been assigned space groups P6 1 22 
(330 atoms per unit cell) and Fddd (704 atoms/cell), 
which have unexpectedly complex crystal structures [5]- 
[T0] . There had been disputes between experimentalists 
and theoreticians regarding the nature of the ground- 
state structure of Mg(BH 4 ) 2 [TTHT4] . To improve the 
reversible hydrogen absorption or desorption kinetics or 
get new metastable polymorphs, a practical method is to 
search for possibility to stabilize the high-pressure phase 
of Mg(BH 4 ) 2 at ambient pressure by advanced theoreti- 
cal or experimental methods. For example, 7-Mg(BH4)2 
is one of the hydrogen-richest solids, reported to be capa- 
ble of storing guest species, such as hydrogen [14]. Most 
recently, new 5, 5' ', and e phases of Mg(BH4)2 were suc- 
cessfully synthesized under pressure [14]. Many of them 
turned out to retain their structure upon decompres- 
sion to ambient conditions. Crystal structures of 7 and 
5 phases were, apparently convincingly, resolved using 
powder x-ray diffraction data obtained with synchrotron 
radiation [14 . Unexpectedly, our phonon calculations 
showed the P4 2 nm structure (S phase) to be dynam- 
ically unstable at ambient pressure, which means that 
the exact crystal structure of 6 phase is still unresolved, 
even for such a simple structure only with 22 atoms per 
cell, and still less for the unknown 8' and e phases [T4] . 
Therefore, the polymorphism and phase diagram of this 
important compound require further investigation. 

Structure searching for ground-state phases of 
Mg(BH4)2 were performed with the ab initio evolution- 
ary algorithm USPEX [T5HT7] . This method has been suc- 
cessfully applied to a wide range of problems [18] [19]. 



When nothing is known about the crystal structure, sep- 
arate searches are usually done with increasing number 
of formula units in the unit cell - as long as computa- 
tional resources allow and until the solution is found. 
We found structure solutions already in searches with 2 
formula units (22 atoms/cell) and 4 formula units (44 
atoms/cell), which yielded structures fully explaining ex- 
perimental results. We did structure prediction at 0, 2, 
5, 10, 15, and 20 GPa. All structures were relaxed by 
the VASP code [20]. During structure searches, we used 
BH4 groups as whole units with fixed bond connectiv- 
ity. All experimentally known structures contain these 
nearly rigid units at pressures below 22 GPa [6 . This 
was done with the Z-matrix representation for the BH 4 
complex anions, keeping the bond connectivity within the 
BH4 tetrahedron, and fully relaxing the bond lengths and 
angles. This approach allowed us to focus searches on 
chemically interesting structures and reduced the search 
space dramatically [17] . By the way, the presence of the 
same nearly rigid BH4 groups, connected to Mg atoms 
by weak bonds implies that zero-point energies of dif- 
ferent structures will be nearly identical and thus can 
be safely neglected when comparing energetics of differ- 
ent structures; our calculations confirm this. We em- 
ployed projector augmented wave method to describe the 
core electrons and the generalized gradient approxima- 
tion of Perdew, Burke, and Ernzernhof (PBE) for ex- 
change and correlation. A cutoff energy of 600 eV, and 
a Monkhorst-Pack Brillouin zone sampling grid with the 
resolution of 2tt x 0.04 A -1 was used. To examine the 
importance of long-range dispersion interactions on this 
compound's thermodynamic stability, the semiempirical 
dispersion-correction method was applied to the ground 
state structures [21] . Phonon dispersion curves were cal- 
culated by using the supercell method as implemented in 
the phonopy code [22] [23] . The powder XRD patterns 
were simulated using the reflex software. 

In our calculations, we found F222 Mg(BEU)2 as the 
most stable phase at ambient pressure, in good agreement 
with previous theoretical results p~2j [13] . The tetrago- 
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FIG. 1. (Color online) (a), (c) Projections of the IAi/acd 
structure along [001] and [010] directions, (b), (d) projections 
of the P4 structure along [001] and [010] directions. Large 
spheres are Mg atoms, small spheres - B (dark or red) and H 
(light or green) atoms. 



nal IA 1 /acd and trigonal P3ml phases are found to be 
the most stable ones in structure searches at 2-5 GPa 
and 10-20 GPa, respectively. Interestingly, within the 
whole pressure range (up to 20 GPa), we did not find the 
PA^nm structure proposed by Filinchuk et al. for the 
S phase [14], but instead found the IAi/acd phase with 
four formula units (44 atoms) per conventional cell and 
Pi phase with 2 formula units per cell at pressures be- 
low 5 GPa (see Fig. 1). Given that the PA^nm structure 
is dynamically unstable at ambient pressure, and based 
on our enthalpy calculations, we hypothesized that the 
IAi/acd and P4 structures might correspond to the ex- 
perimentally observed 5 and S f phases. Further investi- 
gation confirmed this hypothesis, as we will show below. 

Lattice constants of the most interesting candidate 
structures are listed in Table I. The calculated zero- 
pressure lattice constant (a = b = c = 15.79 A) of the 
7 phase is in excellent agreement with the experimen- 
tal value (15.76 A) [14], which gives a benchmark of the 
typical accuracy to expect of DFT simulations for this 



TABLE I. Lattice constants, density (p), bulk modulus (Bo) 
and its pressure derivative (B'), and the total energy of the 
polymorphs of Mg(BH4)2. The energy difference (AE) in- 
cluding van der Waals interactions (AE') is relative to the 
PA^nm structure. Some experimental values (from Ref. 14) 
are also listed for comparison. 



Symmetry 


IAi/acd I4i/amd 


P4 


PA^nm 
(Theory) 


PA^nm 
(Expt) 


a (A) 


7.69 


8.32 


5.58 


5.79 


5.44 


b(k) 


7.69 


8.32 


5.58 


5.79 


5.44 


c(A) 


12.30 


10.52 


5.99 


5.73 


6.15 


p (gem 3 ) 


0.986 


0.984 


0.963 


0.933 


0.987 


Bo (GPa) 


24.0 


31.0 


31.7 


28.5 


28.5 


B' 


4.3 


3.6 


3.6 


3.6 


5.8 


AE (kJ/mol) 
AE' (kJ/mol) 


-17.0 


-15.5 


-13.6 







-21.2 


-19.3 


-15.4 








system. For the IAi/acd structure, the Mg atom occu- 
pies the crystallographic Sb site at (0.5, 0, 0), the B atom 
occupies the 16e site at (0.533, 0.25, 0.375), and the H 
atoms are at the 32g sites with coordinates (0.439, 0.247, 
0.295) and (0.627, 0.123, 0.369). For the P4 structure, 
the Mg atoms occupies the la site at (0, 0, 0) and Id 
site at (0.5, 0.5, 0.5), the B atom occupies at the Ah 
site at (0.25, 0.25, 0.25), and the H atoms are at the Ah 
sites with coordinates (0.459, 0.282, 0.210), (0.139, 0.327, 
0.089), (0.223, 0.038, 0.286), and (0.174, 0.354, 0.416). 
From Table I, one can see that relaxing the experimental 
PA2nm structure, one gets unexpectedly large changes in 
the lattice constants - so large that, in fact, the relaxed 
lattice constants of the P4 structure are the closest match 
to the experimental ones [14 . One has to keep in mind 
that what is called the "experimental" cell parameters 
in many cases is a non-unique result of indexing powder 
XRD spectra, and this is the case here. 

IAi/acd Mg(BH4)2 becomes more stable than the 7 
phase at pressures above 0.7 GPa (Fig. 2). In the room- 
temperature experiment, a pressure-induced structural 
transformation is observed for the porous 7 phase, and 
occurs in two steps: The 7 phase turns into a diffraction- 
amorphous phase at 0.4-0.9 GPa, and then at approx- 
imately 2.1 GPa into the S phase [14 . The calculated 
phase transition pressure from the 7 phase to the pro- 
posed S phase with PA2nm symmetry is 1.2 GPa (the 
corresponding phase transition pressure for P4 phase is 
0.8 GPa), which are in good agreement with the exper- 
imental values (0.4-0.9 GPa). We note a tiny enthalpy 
difference between IAi/acd and P4 structures at pres- 
sures around 1 GPa. As pressure increases to 9.8 GPa, 
the P3ml structure becomes the most stable one, in 
agreement with earlier predictions [6] [10] . Bil et al. [8] 
indicated that it is important to treat long-range dis- 
persion interactions to get the ground state structures 
of magnesium borohydrides correctly. We have exam- 
ined the energetic stability of the considered structures 
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FIG. 2. (Color online) Enthalpy curves (relative to the 7 
phase) of various structures of Mg(BH4) 2 as a function of 
pressure. Enthalpies are given per formula unit. The in- 
set shows the energy per formula unit of P4, I4i/amd, and 
IAi/acd structures including van der Waals interactions at 
zero pressure relative to the P^nm structure. 



through a semiempirical Grimme correction to DFT en- 
ergies, stresses and forces [21 (see the inset of Fig. 2). 
When this correction is included, the IA\/acd and P4 
structures once again come out as more stable than the 
P^nm structure, by 21.2 kJ/mol and 15.4 kJ/mol, re- 
spectively. Energetic stability seems to correlate with 
the degree of disparity of bond lengths and atomic Bader 
charges. The P^nm structure has two inequivalent Mg- 
H distances, 2.26 and 2.07 A, compared to 2.11 and 2.07 
A in the IAi/acd structure, and 2.12 and 2.06 A in the P4 
structure. As we can see, the more homogeneous bond 
lengths, the greater stability. Bader charges show the 
same picture: for H atoms, we find them to be -0.63 and 
-0.59 e in the P4 2 nm structure, -0.63 and 0.62 e in the 
P4 structure, and -0.63 and -0.61 e in the IAi/acd struc- 
ture [24 . More homogeneous Bader charges and bond 
lengths in the IAi/acd and P4 structures correlate with 
their greater thermodynamic stability at ambient pres- 
sure, in agreement with proposed correlations between 
local bonding configurations and energetic stability [T3] . 

Our calculations suggest that the P^nm structure, 
proposed by experiment for the 6 phase, is unstable. This 
implies that either density functional theory calculations 
are inaccurate for this system, or experimental structure 
determination was incorrect. To assess these possibili- 
ties, we simulated the XRD patterns of the IAi/acd and 
P4 structures, and compared them with the experimen- 
tal XRD pattern of the 5 phase at ambient pressure (see 
Fig. 3a). One observes excellent agreement, both for the 
positions and the intensities of the peaks (including both 
strong and weak peaks), of the IA\/acd structure with ex- 
periment [14]. The situation is very peculiar: two struc- 
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FIG. 3. (Color online) Simulated XRD patterns of the 
I4i/acd, P4, and P^nm structures of Mg(BH4)2 with the 
X-ray wavelength of 0.770518 A at ambient pressure (a) and 
0.36814 A at 10 GPa (b) in comparison with the correspond- 
ing experimental results (Ref. 14 and Ref. 6). 



tures, IA\/acd and P^nm, have nearly identical XRD 
patterns, both compatible with the experiment - but one, 
IAi/acd, is the true thermodynamic ground state (global 
minimum of the enthalpy), whereas the other, P^nm, 
is not even a local minimum of the enthalpy (dynami- 
cally unstable structure, incapable of sustaining its own 
phonons). In this situation, the true structure is clearly 
IAi/acd. This example gives a clear real-life example 
of the fact that very different structures can have very 
similar powder XRD patterns, making structure deter- 
mination from powder data dangerous, and in such cases 
input from theory is invaluable. The P4 structure also 
has a rather similar XRD pattern, but the peak posi- 
tions are slightly shifted. Comparison with an indepen- 
dent experimental XRD pattern collected at 10 GPa (Fig. 
3b) shows that the peak positions and intensities of the 
IAi/acd structure are once again in excellent agreement 
with the experimental data [6], while the strong peaks 
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of the P4 structure at 9.9°, 11.6°, and 11.8° obviously 
deviate from the observed ones. This reinforces our con- 
clusion that the IAi/acd structure is the best candidate 
for the high pressure 8 phase. At pressures below 10 GPa 
a mixture of IA\/acd and P4 phases is possible, as the 
XRD peaks of these two structures are quite similar. We 
remind that in the experiment, the 5 and 5' phases are 
nearly indistinguishable [14]. 
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FIG. 4. (Color online) Phonon density of states of (a) the P4 
phase and (b) the IA\/acd phase at ambient pressure. 



Filinchuk et al. demonstrated the bulk modulus of the 
S phase (28.5 GPa) is almost three times higher than 
that (10.2 GPa) reported by George et al. by fitting 
the Murnaghan equation of state [H] . Our third-order 
Birch-Murnaghan [25 fits of the equation of state yielded 
bulk moduli of the IAi/acd and P4 structures equal to 
24 GPa and 31.7 GPa, respectively, consistent with the 
measured value (28.5 GPa) [14]. The observed large den- 
sity difference with respect to the 7 phase at ambient 
conditions (44%) is equally consistent with 45% (43%) 
for IAi/acd (P4) structures [14]. Therefore, it is difficult 
to discriminate between the IA^/acd and P4 structures 
by their compression behavior, density or bulk modulus. 
Our calculations show that the IAi/acd structure does 
not only match all experimental observations for the 5 
phase and has the lowest enthalpy among all sampled 
structures at the relevant pressure range, but is also dy- 
namically stable - phonons were computed at 0, 5 and 
10 GPa. The phonon densities of states (PDOS) of P4 
and IAi/acd phases at ambient pressure are shown in 
Fig. 4, and once again we see a great degree of simi- 
larity. The similarity of all characteristics of these two 
phases parallels the observed similarity of characteristics 
of the S and 5' phases and invites one to propose that 
while the 141/acd structure corresponds to the 5 phase, 
the 5' phase may have the P4 structure. 

In conclusion, we performed a systematic structure 
search for the high pressure phases of magnesium borohy- 
drides and identified two tetragonal structures with space 



groups IAi/acd and P4 as candidates for the unresolved 5 
and 5' phases, respectively, which can be kinetically sta- 
ble upon the decompression to the ambient conditions. 
In particular, the IAi/acd structure is much lower in en- 
thalpy than the earlier reported structures, and the cal- 
culated density, bulk modulus, and the simulated XRD 
patterns are in excellent agreement with the experimen- 
tal results [14]. This example highlights the importance 
of theoretical simulations in establishing crystal struc- 
tures, when only powder XRD data are available: purely 
experimental solution may be dangerous even for simple 
structures, such as the highly symmetric structure of the 
5 phase containing only 6 non-hydrogen atoms in the unit 
cell. 
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